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ABSTRACT 

We present a simple model for electron transport in a possible layer of exotic nuclear clusters 
(in the so called nuclear pasta layer) between the crust and liquid core of a strongly magnetized 
neutron star. The electron transport there can be strongly anisotropic and gyrotropic. The 
anisotropy is produced by different electron effective collision frequencies along and across 
local symmetry axis in domains of exotic ordered nuclear clusters and by complicated effects 
of the magnetic held. We also calculate averaged kinetic coefficients in case local domains are 
freely oriented. Possible applications of the obtained results and open problems are outlined. 
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1 INTRODUCTION 

Neutron stars are thought to consist of the outer crust, inner crust 
and the core (e.g., Shapiro & Teukolsky 1983; Haensel, Potekhin 
& Yakovlev 2007). The outer crust extends to densities p < pnd> 
where ~ (4 - 6) x 10" g cm“^ is the neutron drip density. 
It mostly contains strongly degenerate electrons and atomic nuclei 
(remnants of atoms fully ionized by the huge electron pressure). 
When the density increases, the nuclei become more neutron rich. 

The inner crust (where Pnd < p < p^c, Pcc ~ Po/2 being the 
density at the crust-core interface and po 2.8 x lO" g cm“^ is 
the density of saturated nuclear matter) consists of ultrarelativistic 
highly degenerate electrons, very neutron rich nuclei and free de¬ 
generate neutrons (dripped off the nuclei). With increasing density 
in the inner crust, the fraction of free neutrons (among bound and 
free nucleons) grows up, so that free neutrons give larger contribu¬ 
tion to the pressure, while atomic nuclei become loosely bound as 
if dissolving in the neutron sea. Finally, the nuclei disappear at pcc 
which signals the transition to the liquid core of the star. The crust 
as a whole is about 1 km thick and contains about one percent of 
the neutron star mass. 

The liquid core can extend to about lOpo- The outer core 
(Pcc ^ p 5 2po) consists of strongly degenerate neutrons, with 
some admixture of protons, electrons and muons. The inner core 
(p > 2 Po) may contain hyperons, pion or kaon condensates or free 
quarks or the mixture of these. 

The aim of this paper is to study a possible layer of nuclear 
mantle at the bottom of the inner crust containing various exotic nu¬ 
clear clusters called collectively nuclear pasta. It was predicted by 
Ravenhall, Pethick & Wilson (1983) who investigated properties of 
matter in core-collapsing supemovae. It has been studied by many 
authors as summarized in a discussion below (also see Haensel, 
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Potekhin & Yakovlev 2007 for references to early publications). A 
neutron star mantle can exist or not, depending on specific proper¬ 
ties of nuclear interactions at subnuclear densities. If exist, it forms 
the layer between the ordinary crust and the core, pcm < p < Pcc, 
where p^m ~ 10*"^ g cm“^ is the outer boundary of this layer. The 
mantle is often viewed as a part of the inner crust. Although the 
mantle occurs in a narrow density interval, its mass can be compa¬ 
rable to the mass of the entire crust. As remarked above, the nuclei 
are loosely bound there. In this case nuclear binding becomes as 
low as Coulomb energies of the nuclei. Nuclear pasta forms when 
Coulomb energy is sufficiently strong to rearrange nearly spherical 
atomic nuclei into exotic nuclear clusters. 

The studies of nuclear pasta can be divided into two parts. 
First, one considers hot matter of subnuclear densities for core- 
collapsed supemovae and proto-neutron stars where beta equilib¬ 
rium is absent. The second direction is to investigate colder pasta 
that can be available in neutron stars under the beta equilibrium 
condition. We are interested in the second possibility. Very roughly, 
a transition from the first to the second case occurs at temperatures 
T from about a few times of 10® K to ~ 10'° K. 

There are two major approaches to study nuclear pasta. 

(i) First, one considers one-, two- or three dimensional nuclear 
clusters using different versions of the compressible liquid drop 
model (e.g., Ravenhall, Pethick & Wilson 1983; Hashimoto, Seki 
& Yamada 1984; Lorenz, Ravenhall & Pethick 1993; Watanabe, 
lida & Sato 2000; Watanabe & lida 2003; Nakazato, Oyamatsu & 
Yamada 2009) or Thomas-Fermi theories (e.g., Williams & Koonin 
1985; Oyamatsu 1993; Okamoto et al. 2013; Grill et al. 2014) with 
various mean field interactions. The theory predicts four phases I- 
IV of nuclear pasta in a neutron star mantle. The first phase I, which 
appears with increasing density at p > pcm, is the phase of rods 
(“spaghetti”). Nuclear structures become rod-like there, immersed 
in a liquid of free neutrons. One studies neutron and proton density 
profiles across these rods assuming that their length is large (infi- 
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nite). The phase of rods is followed by phase II of slabs (“lasagna” 
phase). The slabs are filled by nuclear matter and immersed in neu¬ 
tron liquid. One calculates nuclear structure across the slabs while 
the slabs themselves are assumed large (infinite). The phase of 
slabs is followed by phase III of rod-like bubbles (“anti-sphagetti”), 
where free neutrons fill infinitely long bubbles aligned along cer¬ 
tain axis, and the nuclear structures occupy all other space. The 
final phase IV consists of nuclear matter with spherical bubbles of 
neutron (and possibly proton) liquid (“Swiss cheese”). Some cal¬ 
culations predict all four exotic phases, while others predict fewer 
phases or no exotic phases at all. Real rods or slabs are thought to 
be of finite size and occupy some domains. The sizes of nuclear 
clusters and domains as well as orientations of nuclear clusters in 
different domains are uncertain. 

(ii) Second, one simulates the composition and dynamics of 
ensembles of nuclear clusters with various versions of quantum 
molecular dynamics (e.g., Maruyama et al. 1998; Kido et al. 2000; 
Watanabe et al. 2002; Schuetrumpf et al. 2015) or classical molec¬ 
ular dynamics (e.g., Horowitz et al. 2004, 2005; Dorso, Gimenez 
Molinelli & Lopez 2012; Caplan et al. 2014; Horowitz et al. 2015). 
Such simulations predict a variety of possibilities, first of all nu¬ 
clear clusters of one type, e.g. rods, but of finite lengths and differ¬ 
ent sizes. They predict also mixtures of different clusters (e.g., rods 
and slabs). Recently, more exotic structures have been obtained, 
such as slabs with a lattice of bubbles (“nuclear waffles,” Schnei¬ 
der et al. 2014), gyroids (Nakazato, Oyamatsu & Yamada 2009; 
Schuetrumpf et al. 2015), intertwined lasagna and other structures 
which do not resemble usual pasta (Alcain, Gimenez Molinelli 
& Dorso 2014), as well as long lived topological defects (e.g., 
Horowitz et al. 2015). At high temperatures T > 10*® K nuclear 
pasta is mostly in liquid state but at low temperatures (in neutron 
stars) it freezes into ordered or disordered clusters. 

The existence of nuclear pasta does not affect significantly the 
equation of state of the matter and neutron star models. However 
it can strongly influence neutrino emission (e.g., Lorenz, Raven- 
hall & Pethick 1993; Kaminker et al. 1999; Gusakov et al. 2004; 
Horowitz et al. 2005) and transport properties of the neutron star 
mantle (e.g., Pons, Vigano & Rea 2013; Horowitz et al. 2015). We 
will focus on electron transport which is thought to be dominant 
there. To simplify the problem we assume the existence of nuclear 
pasta domains with nuclear clusters of one type in each domain. 
We include the effects of strong magnetic fields which can be very 
important. 


2 ONE DOMAIN 
2.1 General expressions 

We start with the electron transport properties in one domain as¬ 
suming that electron mean free paths are larger than spacing be¬ 
tween nuclear clusters but smaller than domain size and neglecting 
boundary effects. For example, consider ensembles of rod-like or 
slab-like nuclear structures. The electrons are mainly scattered by 
proton charge distributions created by these structures. 

Consider a strongly degenerate and ultra-relativistic electron 
gas which is slightly out of thermodynamic equilibrium due to the 
presence of weak gradients of temperature T and electron chemical 
potential yu as well as due to a weak electric field E. Kinetic prop¬ 
erties of the electrons are calculated from the Boltzmann equation 


for the electron distribution function f(r,p, t). 


df 



( 1 ) 


where v, p and -e are the electron velocity, momentum and electric 
charge, respectively; B is the external magnetic field, and /{/) is 
the collision integral. The distribution function is normalized as 


Ue being the electron number density. Without any loss of generality 
it is sufficient to consider a stationary problem, with dfjdt = 0. 
Because the electrons deviate from equilibrium only slightly, we 
can set f = fo + Sf, where 


/o = 


exp 


kBT) 


+ 1 


(3) 


is the local equilibrium Fermi distribution (e being the electron 
energy and the Boltzmann constant), and 6f is a small non¬ 
equilibrium correction. 

It is well known (e.g., Pitaevskii & Lifshits 1981) that one can 
set / = /o on the left-hand side of (1) in all the terms except for the 
one containing the magnetic force (because the latter term vanishes 
at / = fo). Then 


or op \c / Op 


dfo 

dp 


dp 


(4) 


where I{6f] is the linearized collision integral. 

Since fo depends on r through T(r) and p(r), equation (4) can 
be rewritten as 


vF- 



dp 


im. 


where we have introduced the notations 


(5) 



( 6 ) 


Et being the effective electric field. While calculating the electron 
kinetic coefficients we can treat VT, Vyu, E and E, as arbitrary 
constant small vectors. 

At the next step we should specify the linearized collision in¬ 
tegral I{6f}. In order to study the main properties of kinetic coef¬ 
ficients one can often use the traditional simplified relaxation time 
approximation. 


im = -— = -ySf, (7) 

T 

where t is the relaxation time and v = 1/t is the effective collision 
frequency. In the simplest case when electrons scatter elastically by 
some heavy centers, one has v = n\vao, where u; is the number den¬ 
sity of scatterers and cto is the transport cross section; v, t and 
can depend on the electron energy e. If one interests in a dipole-like 
deviation of the electron distribution function from the equilibrium 
one (which is relevant for studying electric and thermal conductiv¬ 
ities as well as thermopower), it is sufficient to set d/ = v •«, where 
the vector u (to be determined) can depend on the electron energy 
e but not on orientation of its momentum p. Then I[6f} = —vv u. 

The case of nuclear pasta is more complicated because we deal 
with anisotropic medium. As far as nuclear clusters are concerned, 
we assume that we have a symmetry plane (x,y) (perpendicular 
to nuclear rods or plates) and a symmetry axis z (perpendicular to 
the plane). Without any loss of generality we may assume that the 
magnetic field B lies in the (x, z) plane. 
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Table 1. Possible pasta phases in a neutron star mantle 


Number 

Phase 

Leading v 

Weaker v 

I 

Rods 

'"p 

Va 

II 

Slabs 

Va 

'"p 

III 

Rod-like bubbles 

’"p 

Va 

IV 

Spherical bubbles 

Va = Vp 

— 


A natural generalization of the relaxation time approximation 
(7) to this case is 


l[Sf] = -Va - Vp Vp • «p, (8) 

where Va and Vp are generally different effective collision frequen¬ 
cies along the symmetry axis and in the symmetry plane, respec¬ 
tively; Vp = (Va, Vj,); Up = {ux,Uy), and u, as before, can depend 
on e but not on orientation of p. The difference befween Va and 
Vp reflecfs the anisotropy of electron collision frequencies in non- 
spherical nuclear clusters. Such an approximation is often used for 
studying electron transport in anisotropic solids (see, e.g., Askerov 
1981). In the absence of anisotropy, one has Va = Vp and the prob¬ 
lem reduces to a single traditional effective scattering frequency (7) 
in isotropic medium. 

As discussed in Section 1, there could be four phases I-IV of 
nuclear pasta in a neutron star mantle; they are listed in Table 1 
in order of appearance with increasing density at p > pcm. Phase 
I consists of rods. Form symmetry consideration one can expect 
that the most frequent electron momentum transfer occurs when 
electrons move across the symmetry axis (across the rods), so that 
Vp should be larger than Va. We will call Vp the leading collision 
frequency and Va the weaker collision frequency (in phase I). 

Phase I is followed by phase II of slabs (Table 1). Now the 
leading collision frequency Va corresponds to electrons moving 
along the symmetry axis; they can scatter effectively by slabs which 
are perpendicular to the axis. The weaker collision frequency is 
then Vp. 

After slabs one has phase III of rod-like bubbles, where the 
collision frequencies are expected to have the same ordering as in 
phase I, Vp » Va. 

The final phase IV is that of spherical bubbles. In this case 
Va = Vp (electron collisions are isotropic) as in the case of ordinary 
spherical nuclei in the main body of the neutron star crust (at p < 
Pern)- 

For example, according to model I of Oyamatsu (1993), Pmc = 
0.973 X lO*"* g cm“^. The transition to phase II occurs atp = 1.24 x 
lO*'^ g cm“^; to phase III atp = 1.37 x 10*"' g cm“^, and to phase IV 
atp = 1.42x 10*"* gcm“^. Finally, phase IV ends atpcc = 1.43x10*"' 
g cm“^. In a typical neutron star of mass 1.4 Mq and radius about 12 
km the nuclear pasta layer can be about 100 m thick and have mass 
~ 0.01 Mq; it is nearly as massive as the entire overlying ordinary 
crust of spherical nuclei. 

Equation (8) reduces the linearized Boltzmann equation (5) to 

V ■ F - Cl>(v X i>) ■ U = -Va V-U, - Vp Vp • Up, (9) 

where tu = eBcp is the electron gyrofrequency (e includes the 
electron rest-mass energy), and b = BIB is the unit vector in the 
direction of B. In the matrix form 


The solution can be written as 

1 

Ua ^ (11) 


where A = v^Va + a)^(Vpb\ + The matrix Acji, which will be 
denoted as A, is 


A = 


VaVp -I- -w6jVa 

wfojVa VaVp 

<^bj,Vp 


-(jjb^Vp 

Vp -I- (j/bl 


( 12 ) 


where b. = cos 6, b^ = sin 6 and 0 is the angle between B and the 
symmetry axis z. 

This gives us the vector u and f = fo + Sf. Then we can 
calculate the densities of electric and thermal currents, 

j = -T I d^pvf, o = -- f d^p V (e - u) f. (13) 

(27rrr)3 J ' ^ (Intif J 

Formally, j and q have the following structure 


j = &E,+pVT, j = -0TE,-fjVT, (14) 


where hat means 3x3 matrix (tensor) of anisotropic kinetic coeffi¬ 
cients, (T is the electric conductivity tensor, while two other tensors, 
j8 and fj, are auxiliary. The same tensor p enters both expressions, 
for j and q, as a consequence of the Onsager reciprocal relations. 

For applications, it is instructive to rewrite (14) in the form 

E, = 'Rj-a'7T, q = -&Tj-kVT, (15) 

where K = &^\k and a are the tensors of electric resistivity, ther¬ 
mal conductivity and thermopower, respectively. Again, the same 
tensor a enters the expressions for E, and q owing to the Onsager 
relations. The thermal conductivity and thermopower are given by 
k = f] - TP'kp and a = 'R ft, respectively. 

So far our consideration has been quite general, valid for the 
electron gas of any relativity and degeneracy. The analysis is con¬ 
siderably simplified for strongly degenerate electrons. In this case 
the main contribution to electron transport comes from the elec¬ 
trons with energies \e — p\ <k^T in the narrow thermal width of the 
Fermi level. Standard calculations yield 


e^ricA 
cr = -, 


k X fj - 


7p-n,.k^TA 

3ra,A 


( 16 ) 


where A is given by (12), A is defined in (11), m, = pjE is the 
effective electron mass at the Fermi surface (p includes while 

(x>, Va and Vp have to be taken at e = p. Note that a and k satisfy the 
Wiedemann-Franz law. 


3e2 


■ cr. 


(17) 


which is natural for strongly degenerate electrons in the relaxation¬ 
time approximation (e.g. Ziman 1960). 

The tensors & and k have similar structure, determined by the 
matrix A. Therefore, either electric or thermal conduction is defined 
by six different transport coefficients of three types: 


(i) Three diagonal terms (A„, Ayy, A^) are generally different 
reflecting anisotropy of the transport along all three axes; 

(ii) There are two basically different Hall coefficients {A^y = 
—AyP, Ay^ = -Ag.) and two different Hall parameters. 


-*'p 

-aib^ 

0 ' 

/ 

Ux 

' 


' F, ' 

(i>b. 


3 

1 

Uy 


= 

Fy 

0 

cab^ 





1 F, j 


Xp = tu/Vp, (18) 

(10) associated with two collision frequencies, Va and Vp; they describe 
gyrotropic character of electron transport in a magnetic field; 
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(iii) There are two equal coefficients Aj^. = for non¬ 
vanishing fij- and B,; 

(iv) One has A^piB) = A/jai-B), in agreement with Onsager re¬ 
lations. 


Naturally, the electron transport in nuclear pasta is richer in physics 
than in ordinary neutron star matter. 

Let us present also the resistivity tensor 



u)b^ 

0 

-Li)b, 


Cjb: 

0 

3 

1 

Va 


(19) 


Its structure is simpler than that of & or k. 

If & is found from equation (16), one can calculate the ther¬ 
mopower tensor as 


n-k^TK da-ip) 
3e dfj 


( 20 ) 


Such relations are standard (Ziman 1960) for scalar thermopower 
of strongly degenerate electrons in the relaxation time approxima¬ 
tion. Evidently, they remain valid for tensor quantities. In order to 
use them one should express o), v^, and Vp in & as functions of e 
and set e = /r. In addition, one should express through yu and 
differentiate over /j afterwards. 

The presented expressions are sufficient to estimate the elec¬ 
tron transport properties in one domain of nuclear pasta. Below we 
analyze the transport coefficients in limiting cases. Since the struc¬ 
ture of (T and k is the same, we will mainly discuss & meaning that 
one can immediately find k using equation (17). 


2.2 No magnetic field 

In the limit of B —> 0 the transport coefficients are greatly simpli¬ 
fied. 


^ y-l 

P 


0 

0 


( 21 ) 


Here, the conductivity is characterized by the two different kinetic 
coefficients, along and across the symmetry axis. 


CaO - 


v-vn. 


a-pQ 


Vpin, 


( 22 ) 


In phases I and III of nuclear rods or rod-like bubbles one expects 
Vp Va (Table 1). Then the conductivity CTao along the symmetry 
axis is much larger than the leading conductivity CTpo across this 
axis. If we view the leading transport as ‘normal,’ then the trans¬ 
port in the perpendicular direction is abnormally fast. In phase II of 
nuclear slabs one has Va » Vp. The normal leading transport along 
the symmetry axis is much slower than across it (cTao 'K fpo). Fi¬ 
nally, in phase IV of spherical bubbles Va = Vp, so that the transport 
is isotropic (‘normal’ in all directions, as in ordinary neutron star 
crust). 


Then 


m.Vp(l +x^) 


1 

0 


0 

0 




Vp(l + X2)/Va , 


( 23 ) 


The conductivity in the symmetry plane is gyrotropic (governed by 
the magnetic field and the collision frequency Vp), while the con¬ 
ductivity along the symmetry axis is unaffected by B and deter¬ 
mined by Va. 

If B lies in the symmetry plane, along the x axis, then 


m.Vp(VaVp -I- LiV) 


VaVp +ilY 0 0 

0 VaVp -CJVp 

0 tuvp v^ , 


(24) 


The conductivity along B is again unaffected by the magnetic field 
(being naturally determined by Vp), while the conductivity in the yz 
plane (across B) is gyrotropic and regulated by Vp and Va. 

If Va = Vp = Vo, the conductivity is the same as in a basically 
isotropic medium immersed in a magnetic field; see Section 2.6. 


2.4 Three regimes 


If we increase B from 6 = 0, the Hall parameters Xp and Xa grow 
up and the electron transport is affected by the magnetic field. The 
regime of sufficiently small (smaller than 1) Hall parameters corre¬ 
sponds to weakly magnetized electrons, whereas the regime of high 
(higher than 1) Hall parameters to strongly magnetized electrons. 
Because Xp and Xa are not equal, there is also the regime of interme¬ 
diate magnetization in which the electrons are already magnetized 
with respect to one Hall parameter but are still not magnetized with 
respect to the other one. These three magnetization regimes have 
different effects on various conductivity components at different in¬ 
clinations 0 of 6 to the symmetry axis. 

As for the 6-dependence of the conductivity, all conductivity 
components are divided into non-Hall (xx, yy, zz, and zx ones) and 
Hall ones (yx and zy, we do not mention associated non-diagonal 
components like xz or xy). The non-Hall components are even func¬ 
tions of w ex 6 while the Hall components are odd functions. In the 
limit of weak electron magnetization the diagonal components stay 
nearly constant, while cr-^ increases as 6^. All non-vanishing Hall 
components increase proportionally to 6. 

In the limit of high electron magnetization (Xa » 1, Xp » 1) 
we obtain 


m,(vp sim 6 + cos^ 9) 


sin^ 6 -t (XaXp) ' -x^ * cos 6 cos 9 sin 9 


Xj ‘ cos 9 
cos 9 sin 9 


(x„Xp) ' 

xl' sin0 


-Xp ‘ sin 9 
cos^ 9 + Xp^. 


(25) 


With increasing 6 in this regime the non-Hall components either 
stay constant (as cTj-j., ct^ and a-^^) or decrease proportionally to 
6“^, whereas the Hall components decrease as 6“'. The decrease of 
the conductivity components at strong magnetization is evidently 
explained by frequent rotation of electrons around magnetic field 
lines, with many rotations between successive collisions. 


2.3 Finite magnetic field along or across symmetry axis 

If 6 7^ 0, transport properties are generally complicated. The sim¬ 
plest case takes place if B is directed along the symmetry axis. 


2.5 Resistivity 

The structure of the electric resistivity tensor "R, equation (19), is 
simpler than of the conductivity. The diagonal components R-ix = 
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log Xp log JCp log Xp 

Figure 1. (Color online) Components of the electron conductivity tensor (Tap in phases I or III (in units of the leading conductivity CTpo across the symmetry 
axis for = 0) versus the leading Hall parameter Xb = tu/vb oc 5 at three angles 0 = 0, 45 and 90° (from left to right) between B and the symmetry axis. The 
ratio of the effective electron collision frequencies across and along the symmetry axis is taken to be Vp/Va = 10. Shown are all six generally different tensor 
components, cTxx, o’yy, o'zi (thick lines of different types) and (Tyx, and cr^y (thin lines). 





logJta Ios'ka logJCa 

Figure 2. (Color online) Same as in Fig. 1 for electron conductivity tensor in phase II (in units of the leading conductivity Cao along the symmetry axis 
for B = 0) versus the leading Hall parameter Va = cj/Va (proportional to 6) at 6 = 0, 45 and 90° (from left to right) and Va/Vp = 10. 


Kyy = Kj, and = '/?a are totally independent of B, being equal 
to corresponding inverse diagonal components (22) of the electric 
conductivity tensor at S = 0, (Rp = l/Cpo and'??a = l/Cao-The lead¬ 
ing resistivity (iRp or H J is higher than the resistivity in the opposite 
direction. The Hall components = -Kyx and Ky^ - —H^y are di¬ 
rectly proportional to tu oc 6. In the regime of weak magnetization 
the diagonal components are larger than the Hall ones while in the 
the regime of strong magnetization the Hall components are larger. 
Actually, there is only one non-trivial Hall component of K instead 
of two formal Hall components and Ky^ in (19). To prove this 
one can rotate coordinate axes in the xz plane in such a way for the 
z-axis to become parallel to B. After this rotation the tensor would 
contain the only one Hall component = —Ky^ = to¬ 

tally independent of electron collision frequencies Vp and (prov¬ 
ing that the Hall components do not produce any dissipation al¬ 
though produce an important effect of Hall drift of the magnetic 
field). Such a rotation trick does not work out for conductivity ten¬ 
sors - there are generally two non-trivial Hall components of the 
conductivity. 

The simplicity of the electric resistivity tensor is also pro¬ 
nounced in the expression for the Joule heat Q = (j ' ' j) 


[erg s ' cm ^], 

Q = n,jl + njl (26) 

where and j.^ are the components of the electron current j across 
and along the symmetry axis, respectively. 

2.6 Phase of spherical huhhles 

This phase is the simplest one because the medium is basically 
isotropic, Va = Vp = Vq. The conductivity and resistivity tensors are 
especially simple in the coordinate system with the z-axis parallel 
to B. For instance, the electric conductivity and resistivity tensors 
take the well known form 



( 

O-x 

-O-H 

0 




0 ' 

& = 

cth 

cr± 

0 

, R = 

-Rn 

Rx 

0 


. 0 

0 

<^11 , 


. 0 

0 

■^11 , 


where the subscripts || and ± refer to conductivities and resistivities 
parallel and perpendicular to B, respectively, while the subscript H 
labels the Hall terms. The Hall terms in & and 'R have naturally 
























6 D. G. Yakovlev 


different signs. The expressions for these components are 


o-|| 

(Ty 

= 0-0 

1 

1/(1+x2) 


' ^11 ' 

II 

0^1- 

' 1 ' 

1 



_ x/(l +x^) 


, '^H , 

. ^ , 


with CTo = e^n^lvQin,, x = oj/vq being the single Hall parameter in 
the given case. Therefore, iT|| = cto = l/Kw - l/K± are independent 
of B, and Kn oc B. In the regime of weak magnetization, at x <K 1, 
one has nearly constant ^ cr|| but cth ‘x B, while at x » 1 one 
has (Tx oc B^^ and cth ocX b K 


2.7 Illustrative examples 

By way of illustration, in Fig. 1 we show the components of the 
electron conductivity tensor given by equations (16) and ( 12 ), 
in phases of rods or rod-like bubbles versus the leading Hall pa¬ 
rameter Xp, equation (18). The components are expressed in terms 
of the leading (‘normal’) conductivity Cpo at 6 = 0 , equation ( 22 ). 
Recall that the presented ratio cr^^/cr^o is the same as x„^/xpo. The 
three panels correspond to the three inclination angles 0 = 0, 45 
and 90° of B to the symmetry axis. The ratio of the leading and 
weaker electron collision frequencies is taken to be Vp/Va = 10. We 
present all six different tensor components, cr^x, cTyy and cr^ (short- 
dashed, solid, and long-dashed thick lines, respectively), as well as 
o'yx, cr^x (Tg, (dot-dashed, short dot-dash-spaced and long dot- 
dash-spaced thin lines, respectively). 

The conductivities plotted in the left-hand (0 = 0) and right- 
hand (0 = 90°) panels are described in Section 2.3 (equations (23) 
and (24), respectively). At 0 = 0 the conductivity cr^ is unaffected 
by the field, being determined by the weaker collision frequency 
Va. At 0 = 90° the conductivity is not affected by B, being 
determined by Vp. If Xp 0, we reproduce the field-free limit 
(Section 2.2), with the leading (‘normal’) and CTyy conductivi¬ 
ties 10 times smaller than the ‘abnormal’ conductivity cr^ (because 
= 10 , in our example). 

In the limit of weak electron magnetization the diagonal com¬ 
ponents stay nearly constant, while grows as B- and all non¬ 
vanishing Hall components increase as B. In the limit of high elec¬ 
tron magnetization the non-Hall components either stay constant 
(as (Txx, cTzz and cr-j at 0 = 45° in Fig. 1) or decrease as B^^ (like 
(Tyy in the same case). In this regime the Hall components decrease 
as B '. 

Fig. 2 plots the same conductivities as in Fig. 1 for 0 = 0, 
45 and 90° but in phase II of nuclear slabs. The conductivities are 
expressed in units of the leading (‘normal’) conductivity cr^o at B = 
0, equation (22); they are shown as functions of the leading Hall 
parameter x^. The ratio of the leading to weaker electron collision 
frequency is again assumed to be Va/Vp = 10. The main features 
of the conductivities are similar to those in Fig. 1 for the phase of 
rods. In the regime of weak magnetization for slabs we also have 
the abnormal conductivities (cTa-x and cr^) much higher than the 
‘normal’ conductivity (cr^). 


3 AVERAGING OVER DOMAINS 

In Section 2 we have considered the electron transport in one do¬ 
main of exotic nuclear structures. This transport is described by 
3x3 matrices of kinetic coefficients; the matrix elements be¬ 
have differently and strongly depend on the magnetic field B and 
anisotropic properties of nuclear structures. Macroscopic equations 


of electron heat and charge transport contain the conductivity (re¬ 
sistivity) tensors averaged over an ensemble of domains. In prin¬ 
ciple, the averaging should take into account correlation properties 
of electron transport in different domains. However, the theory of 
nuclear pasta phases is still not very elaborated (Section 1), so that 
no reliable information has been obtained about predominant ori¬ 
entations and other statistical properties of the domains. 

Therefore, it seems instructive to consider the simplest model 
of freely oriented domains. We take the conductivity tensors ob¬ 
tained in Section 2 and average them over free orientations of the 
domains. It is evident, that after such averaging we obtain the con¬ 
ductivity tensors for basically isotropic medium in a magnetic field. 
In a coordinate frame with the z-axis along B, the conductivity or 
resistivity tensors should have the same form (27) as in the phase 
of spherical bubbles (Section 2.6). Any tensor is characterized by 
three components which describe transport along and across B as 
well as the Hall transport, although these components are generally 
more complicated than (28). Therefore, it is sufficient to determine 
these three components. 

To this aim, let us take, for instance, the electric conductivity 
tensor (16) and calculate cth = (b ■ & ■ b), cry = {e\ ■ & ■ e\) and 
Ch - € 2 ), for a given orientation of the symmetry axis. 

Here, ei and 62 are two orthogonal unit vectors perpendicular to b, 
with €2 = €[ X b. We obtain 


/ '1 

o-\\ 

^2 

VaVp/jJ -1- + (1? 

O-y 

e lie 

VaVpCel + ep + 

, O-H , 


, VaWiJ -1- Vpiobj + ei-e2c(v^ - VpVa) , 


Now the averaging over orientations of nuclear clusters, which we 
denote as (...), is equivalent to averaging over rotation of the (ei, 
62 ) plane about the b axis and to subsequent averaging over orienta¬ 
tions of the b axis. These averagings can be done analytically with 
the result that in any reference frame with the z axis parallel to B 


(&) = 


(o-x) -<o-H> 0 

((fn) (o-y) 0 

0 0 <cr||) 


The averaged conductivity components are given by 


( <^ll> 1 


( {(1? + r^)(w^ + VpVa)H - ''pi 

(0-x) 

m,cj^ 

{ (vaVpftU^ - \r^)H + Vp) 

1 <0-H> J 


( w - (x)v.y-H 


where 


H = 


(js) 'arctan(j/r) 

( 2 rsy^ In [(r -t- i)/(r - x)] 
l/(v’ + (xl^Va) 


at Va > Vp, 
at Va < Vp, 
at Va = Vp, 


(30) 


(31) 


(32) 


with r = + VaVp) and s - o) -^|vp - Va|. 

All three conductivities (cr||), {(Ty), and (cth) depend generally 
on B (or on ai) and on collision frequencies Va and Vp. In the limit 
of weak electron magnetization we naturally have (cr||) = {(Ty) = 
{(To), where 


{ 0 - 0 ) = 


e'«e<v-‘> 


<v-‘) = 


3v„ 3Va 


(33) 


(cTo) being the field-free averaged conductivity, and (v“*) the aver¬ 
aged inverse collision frequency (averaged relaxation time) of elec¬ 
trons; {(To) is given by the familiar Drude formula. According to 
(29) and (31), in the same limit of weak electron magnetization the 
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Figure 3. (Color online) Electron conductivities (ct/) (averaged over orientations of nuclear clusters, in units of the leading conductivities (22) at B = 0) along 
B {i =11, solid lines), across B {i =_L, short-dashed lines), as well as the Hall conductivity (/ =H, long-dashed lines) as functions of the leading Hall parameter 
for the phases of rods, slabs and spherical bubbles (from left to right). For rods and slabs the conductivities are given for two ratios of the leading to weaker 
collision frequencies, 10 and 100 (thick and thin lines, respectively). See text for details. 


Hall term behaves as 


(o-h) ~ 


e^n^iv) 0) 

-T- OC B, 


(34) 


where (v) = (Va + 2vp)/3 is the angle averaged collision frequency. 

For illustration, in Fig. 3 we show all three electron conduc¬ 
tivities (cr,) averaged over orientations of domains, in units of the 
leading (‘normal’) conductivities at 6 = 0. The conductivities 
along B (i =||) are plotted by solid lines, across B (i =±) by short- 
dashed lines, and the Hall conductivities (i =H) by long-dashed 
lines. These conductivities are presented as functions of the lead¬ 
ing Hall parameters (18) for the phases of rods, slabs and spherical 
bubbles (left, middle and right panels, respectively). For the phases 
of rods and slabs the conductivities are presented at two ratios of 
the leading to weaker collision frequencies, 10 and 100 (thick and 
thin lines, respectively). 

As discussed above, in the weak magnetization regime the 
conductivities (cr||) and (ctx) are almost equal and independent of 
B. At large difference between v.^ and Vp these conductivities are 
much larger than the corresponding ‘normal’ conductivities Cao or 
(Tpo (Section 2, equation (22)). This results from the presence of 
abnormally high conductivities in any domain. By averaging over 
orientations of domains we inevitably allow for abnormally high 
transport in any direction. 

In the regime of strong magnetization the longitudinal con¬ 
ductivity, again, becomes independent of B. However, it depend on 
Va/Vp as 



(35) 


where 


S = 


(yC) ‘arctan(f/ 7 ) 

(2yO-' ln[(7-H^)/(r-^)] 

1 /Va 


at Va > Vp, 
at Va < Vp, 
at Va = Vp, 


(36) 


with 7 = yv^ and f = y|vp - Va|. The higher the difference be¬ 
tween Va and Vp, the lower (cr||) (measured in units of the leading 
field-free conductivity). 

In phase of spheres (cr||) is naturally independent of B being 
equal to (cto) at any B. 


As for the transverse averaged conductivity at high magneti¬ 
zation, it is suppressed with growing B as (cTj^) oc B^^ (Fig. 3). The 
appropriate asymptote is 


{o-±) = 


e^UcVp 

2 m,oj~ 


(1 + VaS). 


(37) 


The averaged Hall conductivity at high magnetization behaves 
as (cth) « It is given by the universal asymptote 

2 

e He 

(o-h) ^ , (38) 

m,u) 

which is independent of collision frequencies (meaning that (cth) 
becomes collisionless) and of specific pasta phase. In this regime, 
(cth) is governed by classical Hall drift of electrons in crossed elec¬ 
tric and magnetic fields (e.g., Pitaevskii & Lifshits 1981). Gener¬ 
ally, the anisotropy of collision frequencies has larger effect on av¬ 
eraged conductivity components in phase of slabs than in phase of 
rods. 

Finally, let us discuss the tensor of averaged resistivity {K) 
which has the same form as (27). A simple averaging of the resis¬ 
tivity tensor (19) gives 


(^||> = (7?r) = (%> = (^H> = 


(39) 


This tensor has the same structure as the tensor (28) in phase of 
spheres. The only collision frequency involved here is the averaged 
collision frequency (v). It determines the averaged longitudinal and 
transverse resistivities which are equal and independent of B. The 
Hall resistivity C/^h) exceeds (7?o) at high magnetization, tu > (v). 

Let us stress, that, generally, (fk) ^ (<t)“‘ and the difference 
can be substantial in anisotropic media. While calculating Joule 
heat (26) one should definitely use the averaged resistivity, but not 
the inverse averaged conductivity. Large difference between Va and 
Vp does not produce any significant difference between (;Ro) and the 
leading non-averaged resistance (l/iTao or l/cTpo). 


4 COLLISION FREQUENCIES 

So far, we have treated the collision frequencies Va and Vp as given. 
Let us discuss how to calculate them. 
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Figure 4. (Color online) Electric and thermal conductivities (thick solid and 
dashed lines, respectively) versus temperature in phase IV of spherical bub¬ 
bles at p = 1.423 X lO'^ g cm“^. Thin solid and dashed lines show the 
effective collision frequencies v^- and v, for electric and thermal conductiv¬ 
ities, respectively. See text for details. 


Figure 5. (Color online) Reduced dimensionless thermopower components 
ay, a± and tfu (solid, long-dashed and short-dashed lines, respectively) 
given by (40) versus the Hall parameter x = ljIvo- in phase IV of spher¬ 
ical bubbles at p = 1.423 X 10*'* g cm“^ and T = 3 X 10* K. See text for 
details. 


4.1 Spherical bubbles 

The formalism for calculating the collision frequency vo (= Va = 
Vp) in phase IV of spherical huhhles is the same as that for spherical 
nuclei in the ordinary neutron star crust. For illustration, consider 
model I of Oyamatsu (1993), according to which spherical buh- 
bles exist in the narrow density range (from about 1.42 x 10*'* to 
1.43 X 10*'* g cm“*. Section 2.1), just near the interface between the 
mantle and neutron star core. Analytic fits for proton and neutron 
density profiles in nuclear bubbles are given by Haensel, Potekhin 
& Yakovlev (2007) in their Appendix B. 

For certainty, we take p = 1.423x 10*'* g cm“* (i.e., the baryon 
number density «(, = 0.08517 fm^**) although the results are almost 
independent of density within such a narrow density range. The ra¬ 
dius of a spherical Wigner-Seitz cell around one spherical bubble 
is 7?vvs = 15.4 fm. The cell contains ~ 1300 nucleons includ¬ 
ing Ztot ~ 45 protons. They are mainly free neutrons distributed 
almost uniformly. The neutron and proton density profiles within 
the cell decrease near the cell’s center. This creates neutron and 
proton holes which can be collectively treated as a nuclear hole. 
The minimum neutron density in the center is only slightly lower 
than the average neutron density, while the central proton density 
drops to zero. The proton hole is equivalent to vacancies of Zh ~ 11 
protons; whereas the entire nuclear hole is equivalent to vacancies 
of A„uc ~ 48 nucleons. The root-mean-squared radius of the proton 
hole is about 8.2 fm; the neutron hole is shallow but larger. These 
nuclear holes behave as charged quasiparticles and form a strongly 
coupled Coulomb system. The electrons are distributed almost uni¬ 
formly and scatter off charge fluctuations of proton vacancies. The 
electron chemical potential is p 87.3 MeV. The electron trans¬ 
port properties can be calculated by applying standard formalism 
of electron transport in a plasma of spherical nuclei (Potekhin et al. 
1999; Gnedin, Yakovlev & Potekhin 2001; Potekhin et al. 2013) 
and replacing the nuclei by the nuclear holes. 

Fig. 4 shows the temperature dependence of the electric and 


thermal conductivities, a and k (thick lines), in the phase of spheri¬ 
cal bubbles atp = 1.423x10*'* g cm“*. The conductivities are calcu¬ 
lated with the method of quasi-potential for electron-nucleus scat¬ 
tering suggested by A. Y. Potekhin (Potekhin et al. 1999; Gnedin, 
Yakovlev & Potekhin 2001) and modified to consider electron scat¬ 
tering by nuclear holes. The calculations include Coulomb coupling 
of nuclear holes, Debye-Waller factor, multi-phonon processes, 
and shape of spherical bubbles. The results are more advanced than 
those based on the simplified model (7) of the collision integral 
in the relaxation time approximation. The displayed conductivi¬ 
ties are solely produced by electron-hole scattering; other scatter¬ 
ing mechanisms (e.g., the electron-electron scattering, Shternin & 
Yakovlev 2006) are disregarded. Low-temperature suppression of 
the electron-phonon scattering due to electron band structure ef¬ 
fects (Chugunov 2012) is not included. 

The conductivities a and k in Fig. 4 are expressed in CGS 
units (in s“' and erg cm“* s“‘ K“*, respectively). They are remark¬ 
ably close to those calculated for standard models of inner neutron 
star crust of spherical nuclei (without nuclear pasta, e.g. Gnedin, 
Yakovlev & Potekhin 2001; Potekhin et al. 2013; note that the con¬ 
ductivities cr in fig. 7 of the latter paper are plotted - from top to 
bottom - for log T =8, 8.5 and 9, respectively). Two thin lines in Fig. 
4 present the effective collision frequencies Vo- and (expressed in 
s“*) which describe the electric and thermal conduction, respec¬ 
tively (see, e.g., Potekhin et al. 1999). At T > 10® K these fre¬ 
quencies almost coincide. This means that the collision integral in 
the Boltzmann transport equation can be accurately described in the 
relaxation time approximation, equation (7). At lower temperatures 
Vn becomes higher than and the relaxation time approximation is 
less accurate. However, at temperatures T > 10* K, most important 
for many applications, this approximation can still be used, as least 
for semi-quantitative analysis. 

Having the effective collision frequency, we can easily esti¬ 
mate magnetic fields 6 = Sh which correspond to the Hall param- 
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eter x = (jjjv = 1. They are characteristic magnetic fields which 
affect the electron transport. For T = 10* and 10® K with v = 
we obtain Sh ~ 2 x 10‘* and 2 x 10** K, respectively. The B fields 
expected in neutron star interiors can be much higher. The lower 
the temperature, the lower the collision frequency, and the lower 
the fields which affect the electron transport. 

So far we have not discussed the thermopower tensor a. For 
nuclear bubbles, we can present this tensor in the same form as &, 
equation (27), and write a = [n-k^Tl{'iep^vpj\'a, where a is the 
reduced (dimensionless) thermopower and vp ~ c is the electron 
Fermi velocity. From equation (20) one obtains (Urpin & Yakovlev 
1980) 

P xV 

a\\ = 'i+P, ttj. = 3 + --ttH = 1-j, (40) 

1 + X* 1 + X* 

where 'P = 51n(evo-(e))“*/31np|p=pp and, again, x = tu/v^p = 
B/Bh- Note that the expression for ay (which is independent of 
B) is valid as long as the electrons are strongly degenerate, regard¬ 
less of the applicability of the relaxation time approximation (7). In 
contrast, the expressions for aj_ and depend on B through x and 
require the approximation (7). 

Fig. 5 presents all thermopower components for the phase of 
nuclear bubbles at p = 1.423 x 10*"* g cm“* and T = 3 x 10* K as a 
function of x = BIBu, with Ba x 2.03 x 10** G. The thermopower 
ay along the magnetic field is shown by the solid line, while 
and f?H are plotted by the long-dashed and short-dashed lines, re¬ 
spectively. Recall that the thermopower components, contrary to 
the components of &, k and 'R, can change sign depending on elec¬ 
tron scattering mechanism. Signs of the thermopower components 
are regulated by the parameter P. As shown by Urpin & Yakovlev 
(1980), neglecting the effects of the Debye-Waller factor and finite 
sizes of atomic nuclei, one has y’ = -1 for scattering of ultrarela- 
tivistic degenerate electrons in high-temperature crystals (when the 
relaxation time approximation is strictly valid). Indeed, we repro¬ 
duce P = as. T = 10® K when the relaxation time approxima¬ 
tion is really accurate (see above). However, taking into account 
the Debye-Waller factor and finite sizes of the nuclei, we obtain 
P X 1.1 (instead of -1). Therefore, at the bottom of the neutron 
star crust the effects of the Debye-Waller factor and of the shape 
of proton (proton hole) charge distribution are very important. If 
we decrease the temperature from 10® to 10* K, the relaxation time 
approximation becomes less accurate and P increases slowly from 

1.1 to 1.3. 

According to Fig. 5, the magnetic field modifies but not 
very strongly. However, it significantly affects an which regulates 
the Hall drift of the magnetic field. At x < 1 the Hall thermopower 
directly determines the drift velocity (Urpin & Yakovlev 1980). Ac¬ 
tually, the thermopower can be large enough to affect generation 
and evolution of magnetic fields in hot neutron stars. 

4.2 Non-spherical nuclear structures 

Calculations of Va and Vp for phases I-III of non-spherical nuclear 
structures are much more complicated. In principle, they can be 
performed using the technique similar to that for spherical nuclei. 
To this aim one needs detailed studies of proton charge fluctua¬ 
tions and their effects on electron scattering. For instance, useful 
information on vibration properties of nuclear clusters can be ob¬ 
tained using their elastic constants (e.g., Pethick & Potekhin 1998). 
Proton charge fluctuations for the phase of slabs have been ana¬ 
lyzed briefly by Watanabe, lida & Sato (2000). Classical molecu¬ 
lar dynamics simulations of dynamical behavior of nuclear pasta 


have been conducted by Horowitz et al. (2005). However, these 
results seem insufficient for a reliable determination of Va and Vp. 
The required physical ingredients are the effective structure factor 
S(q), Debye-Waller factor and nuclear form-factor F(q). Because 
nuclear clusters are non-spherical, these quantities will strongly de¬ 
pend on orientation of wave vector q, leading to essential difference 
between Va and Vp. To the best of our knowledge, the collision fre¬ 
quencies Va and Vp for non-spherical clusters have not been calcu¬ 
lated. 

In this situation, it seems reasonable to assume that the leading 
effective frequency, Va or Vp, is close to that for models of neutron 
star matter with spherical nuclei, while the other (abnormal) effec¬ 
tive frequency is smaller. The ratio of the leading to the smaller 
frequency can be considered as a free parameter, ranging within 
reasonable limits, for instance, from 10 to 100. It would be difficult 
to be more exact using existing calculations. For illustration, one 
may look at fig 3 of Shternin & Yakovlev (2006) which shows par¬ 
tial thermal conductivities of electrons at 6 = 0 in a neutron star 
crust of spherical nuclei. Assume the presence of nuclear pasta at 
the crust bottom and take the leading collision frequency the same 
as for electron-nucleus scattering in fig 3. Let the weaker collision 
frequency be limited by electron-electron collisions (whose con¬ 
tribution is also plotted in fig 3). Then in the temperature range 
10* - 10® K the ratio of the leading to weaker collision frequen¬ 
cies would be ~ 10. Accurate microscopic calculations of and Vp 
would be highly desirable. 


5 CONCLUSIONS 

We have analyzed electron transport in possible nuclear pasta 
phases in a mantle of magnetized neutron star. Specifically, we have 
performed model calculations of the tensors of electric and thermal 
conductivity, & and k, thermopower, a, and electric resistivity R. 
The model assumes the presence of domains containing (quasi)- 
ordered exotic nuclear structures of one type (rods, slabs, rod-like 
bubbles, or spherical bubbles; Table 1). The collision integral in 
the electron transport equation is described in the generalized relax¬ 
ation time approximation (8) with two, generally different, effective 
collision frequencies, Va and Vp, along and across the symmetry 
axis of corresponding nuclear clusters. For non-spherical nuclear 
clusters one can expect the presence of the leading (normal) col¬ 
lision frequency (va or Vp) and abnormally weak frequency (Table 
1), whereas for spherical bubbles Va = Vp. Our model is simple but 
it allows us to include the effects of strong magnetic fields which 
would be difficult to do using more complicated models. 

In Section 2 we have studied the electron transport in one do¬ 
main and have found it rather complicated, with six generally dif¬ 
ferent components of conductivity tensors. 

In Section 3 we have considered electron transport in matter 
containing freely oriented domains of one type. The tensors of ki¬ 
netic coefficients (30), (31), and (39) averaged over orientations 
of the domains are calculated analytically and have much simpler 
structure than their non-averaged counterparts. Each tensor con¬ 
tains only three different components, along and across B, as well 
as the Hall component. These components show sufficiently sim¬ 
ple dependence on B. Nevertheless, this dependence is less trivial 
than in ordinary basically isotropic media with Va = Vp. For in¬ 
stance the angle-averaged longitudinal conductivities (iT||) and (^cn) 
depend now on B and can be much higher than the leading normal 
conductivities (22) for 6 = 0 determined by the leading collision 
frequencies. This enhancement of averaged longitudinal conduc- 
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tion is produced by the contribution of abnormally high conduction 
with lower collision frequency. The enhancement contradicts some 
expectations (e.g.. Pons, Vigano & Rea 2013; Horowitz et al. 2015) 
that nuclear pasta suppresses the conduction. It seems both possi¬ 
bilities (suppression and enhancement) can be realized; see below. 

In Section 4 we have outlined how to calculate the collision 
frequencies Va and Vp. The collision frequencies and transport prop¬ 
erties in phase of nuclear bubbles can be calculated rather reliably; 
they appear similar to those in traditional phase of spherical nu¬ 
clei. However, the collision frequencies for non-spherical nuclear 
clusters are uncertain and require further attention. 

Much work should be done to clarify transport properties of 
nuclear pasta. The major problem is to reliably calculate effective 
collision frequencies Va and Vp for non-spherical nuclear clusters 
I-III (as outlined in Section 4.2). The problem requires accurate 
studies of proton charge fluctuations in these structures. 

We have treated nuclear clusters in one domain as ordered and 
having the same shapes. If the ordering is not exact and the param¬ 
eters (e.g., Va and Vp) strongly vary, an additional averaging over 
disorder and parameter distributions has to be performed. We have 
assumed that orientations of domains are free. If not, more compli¬ 
cated averaging over orientations is needed. 

We have neglected quantum oscillations of electron kinetic 
coefficients in very strong magnetic fields due to population of 
new electron Landau levels with increasing density (e.g., Potekhin 
1999). If the electrons occupy many Landau levels, these os¬ 
cillations are not too strong and our results should reproduce 
oscillation-averaged kinetic properties. In addition, since nuclear 
clusters are loosely bound, very strong magnetic fields may influ¬ 
ence proton charge fluctuations and affect transport properties in 
this way. All these problems are beyond the scope of this paper. 

Another, more serious problem would be to account for pos¬ 
sible impurities and defects in nuclear pasta. It seems that these 
effects can greatly enhance or suppress electron transport. 

Finally, let us stress, that there are other heat and charge trans¬ 
port mechanisms in nuclear pasta which should be taken into ac¬ 
count along with the mechanisms studied above. For instance, one 
should bear in mind thermal conductivity due to electron-electron 
collisions (e.g. Shternin & Yakovlev 2006) and due to superfluid 
phonons (Aguilera et al. 2009) associated with the presence of free 
superfluid neutrons. 

The results of the present studies can be useful for model¬ 
ing of thermal relaxation between crust and core in young neutron 
stars and also in accreting neutron stars in X-ray transients with the 
crust overheated during accretion episodes (e.g.. Pons, Vigano & 
Rea 2013). In addition, the results can be used to study evolution 
of magnetic fields in neutron stars and related phenomena (e.g., 
Horowitz et al. 2015). 
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